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1 Introduction 



The understanding of the structure of the perturbative series has advanced 
considerably over the recent years, see the review [1]. It became clear that 
perturbative series are at best asymptotic, not even Borel-summable. Based 
on an analysis of singularities in the Borel plane, one can obtain the behaviour 
of the perturbative series for large L, where L is the order of perturbation 
theory. The nearest singularity determines the leading asymptotic behaviour. 

Most of the previous investigations use the large-/3o limit, whose relation to 
the real QCD is unclear. At the first order in l/jSo, singularities in the Borel 
plane are simple poles. At the higher orders, they become branching points. 
However, there is an approach [2,3] based on the renormalization group, which 
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yields results with a real basis in QCD. Singularities in the Borel plane are 
branching points, whose powers are determined by the relevant anomalous 
dimensions, but normalization factors cannot be calculated. 

Effective field theories make use of the fact that a large scale is present, and 
physical quantities can be expanded in inverse powers of this large scale. In 
Heavy Quark Effective Theory (HQET, see the textbook [4]), this scale is the 
heavy quark mass m. Renormalon singularities in HQET were investigated 
in [5,6]. Unhke in QCD, the HQET heavy-quark self-energy has an UV renor- 
malon at positive u, namely u — ^, which leads to an ambiguity in the residual 
mass term. 

A typical matrix element in the full theory, QCD, is expanded in 1/m: 

<^>='^<>^^j:^'<0'>^o{^) (1.1) 

(see (2.8)), with short-distance matching coefficients C, Bj,. . . and long-distance 
HQET matrix elements <J>, <0i>,. . . 

The QCD matrix element <j> contains no renormalon ambiguities, if the op- 
erator j has the lowest dimensionality in its channel. ^ In HQET, we separate 
short- and long-distance contributions. In schemes without strict separation 
of large and small momenta, such as MS, this procedure artificially intro- 
duces infrared renormalon ambiguities in matching coefficients and ultraviolet 
renormalon ambiguities in HQET matrix elements. When calculating match- 
ing coefficients C,. . . , we integrate over all loop momenta, including small 
ones. Therefore, they contain, in addition to the main short-distance contri- 
butions, also contributions from large distances, where the perturbation the- 
ory is ill-defined. They produce infrared renormalon singularities, factorially 
growing contributions to coefficients of the perturbative series, which lead to 
ambiguities ~ (Aqcd/'^ti)" in the matching coefficients C,. . . Similarly, HQET 
matrix elements of higher-dimensional operators <Oi>,... contain, in addi- 
tion to the main large-distance contributions, also contributions from short 
distances, which produce ultraviolet-renormalon singularities. They lead to 
ambiguities of the order Aq^^j^ times lower-dimensional matrix elements (e.g., 
<J>). These two kinds of renormalon ambiguities have to cancel in physical 
full QCD matrix elements <j> (1.1) [7,8]. 

Although this has been shown explicitly only in the large-/9o limit, it is as- 
sumed to hold beyond this approximation. Based on this assumption, one 
may obtain additional information on the structure of the infrared renor- 
malon singularities of matching coefficients, based on ultraviolet renormalons 

^ Otherwise, there may be several ultraviolet renormalons on the positive half-axis, 
leading to ambiguities of the order Aqq^ times lower-dimensional matrix elements. 
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in higher-dimensional matrix elements, which are controlled by the renormal- 
ization group. This modcl-indcpcndcnt approach was applied to some simple 
HQET problems: the heavy-quark pole mass [9] and the chromomagnetic- 
interaction coefficient [10]. 

In the present paper, we investigate heavy-light quark currents. The asymp- 
totic behaviour of the pcrturbative scries for the leading QCD/HQET match- 
ing coefficients (due to the nearest infrared renormalon) was studied in [7,11,12] ^ 
in the large-/9o limit. Here we go beyond this approximation, by using the 
renormalization-group based method. The results on 1/m expansions of QCD 
heavy- light currents arc collected in Sect. 2. We show that the asymptotic 
behaviour of the pcrturbative series for the matching coefficients for all cur- 
rents follows from just four distinct cases, two spin-0 currents and two spin-1 
ones. These cases are considered in Sects. 3 and 4 in detail. Ratios of meson 
matrix elements, such as are given by the ratios of the correspond- 

ing matching coefficients at the leading order in 1/m. The asymptotics of the 
pcrturbative series for this ratio is discussed in Sect. 5. The large two-loop cor- 
rection in this ratio was observed in [11]; here we present model-independent 
results for higher orders which continue this trend. 



2 Heavy-light currents in QCD and HQET 

Heavy Quark Effective Theory (HQET, see the textbook [4]) has greatly ad- 
vanced our understanding of many problems in heavy quark physics. Its La- 
grangian is [13,14,15] 

L = Kiv -DK^^ [Ok + CMOM] + C»(l/m2) , (2.1) 
Im 

Ok — —hyD']_hy , = -hvGai30'°'^hy , 

where hy — fh^ is the heavy-quark field, and Dj^ — D — v{v ■ D). Due to 
reparametrization invariance [16] , the kinetic-energy operator Ok is not renor- 
malized, and its coefficient is unity to all orders in perturbation theory. 

The chromomagnetic- interaction coefficient Cm(/i) can only be computed per- 
turbatively, by matching the amplitudes of an appropriate scattering process 



^ Note a typo in (4.8) of [11]: denominators of both terms with a should be 27r, not 

TT. 
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in QCD and HQET. Solving the renormalization-group equation, we obtain 
CM^C^i^^j K_,Ja,{n)) , (2.2) 



where as{iJ,) is the QCD couphng with ni hght flavours in the MS scheme, 



and for any anomalous dimension ^{as) = 7oQ!s/(47r) + 7i(Q!s/(47r))^ + • • • we 
define 

f \ ?Y 7("s) lo\da, T . 7o (li Pi\ as . .^oA 



The anomalous dimension of the chromomagnetic operator Om is [14,15,17,18] 
= 2Ca^ + ^Ca{17Ca - ISTpni) (^^ j + • • • (2.4) 



The full one-loop correction to has been calculated in [14], and the two- 
loop correction in [18]. It is convenient to choose 

Ho = e^^/^m. (2.5) 
We have (see [10]) 



Cm^l + Cml^^ + ---: (2.6) 

cmi = 2Cf + ^Ca - (sCf + ^Ca) ^ + (IIC^ + 7Ca) ^ • 



Operators of full QCD are expanded in 1/m; coefficients of such expansions 
are HQET operators with appropriate quantum numbers. In the present paper 
we shall consider heavy-light quark currents jo = Qo^Qo = '^f('^s(A*'))i(A*')' 
where F is an antisymmetrized product of n Dirac 7 matrices 

r = 7t^---7r' or 7f •••7r-^V, (2-7) 
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which commutes or anticommutes with f: 



fV = aF^ , cr = ±1 , 



and 7^ = 7" — fjv'^. Here a'^fJ') is the QCD couphng with Uf = n/ + 1 
flavours. At the leading order in 1/m, we have only a single HQET current 
Jo = Qo^hyQ = Z{as{n))j{fj.), while to subleading order in 1/m we get 



(2.8) 



The solution of the renormalization- group equation for Cr{fJ>',ij) is 



Cr{n',f^) = Cr ( #M r ( ^ ) K'da',{^^'))K_,{aM) ■ (2-9) 



«s()"o). 



Here K' involves the n/-flavour /3-function The anomalous dimensions of 
J and J are 



(2.10) 



/So 5\ ^ /2 2 49\ ^ 10^ 1 



7[^ = -2(n-l)(n-3)C^^ 



X U + 



i(5(n - 2)2 - 19)Cf - |(3(n - 2)^ - 19)Ca 



An 



-i(n-l)(n-15)C^/5^(^)' + 



(2.11) 



where n is the number of 7 matrices in F (2.7) (of course, the anomalous 
dimension of the vector current vanishes to all orders). The anomalous di- 
mension 7 [19,20,21,22] of the HQET current j does not depend on F (the 
three- loop term can be found using the method of [23]). The anomalous di- 
mension 7p of the QCD current is known to two [11] and three loops [24]. The 
full one-loop corrections to Cr{m, m) were obtained in [13], and two-loop ones 
in [11,25]. The renormalization-group invariants Cr are given by perturbative 
series in asdio): 



Cr 



L=l 



An 



(2.12) 
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+ 
+ 



.4^2 23 
3" ^4 



+ + 8) C^] - - I (IIC^ + ICa) 
I {{n - 2f - 5) - I (3(n - 2)^ - 4) Ca] {n -l){n-3)j^ 



+ {llCp + 7Ca) {n - l)(n - 3)- . , 



PI? 



where r] = — 1)". It will be possible to find C2 from the known two-loop 
results for Cr{m,m) when 72 will be known. In this paper, our aim is to 
investigate the behaviour of the coefficients at L » 1. 



There arc various prescriptions for handling 75 in dimensional regularization. 
Multiplying T by the anticommuting 7^'-' does not change Cp. This is not true 
for the 't Hooft-Veltman 7^^^. QCD currents with 7^^ and 7^^^ are related 
by finite renormalization factors [26,27] : 

(«T7-0)^, = A-;,^„-,;^„, (aU/)) («T7r<3),, . (2.13) 



where the anomalous dimensions 7p^Ac = 7f ^r^^^ differ starting from 
two loops. In HQET, both currents have the same anomalous dimension 
7, and hence the corresponding renormalization factor is unity. Therefore, 
C-p^nv {fj,' , fj,) differs from Cr(/u',/i) only by K'^^{a'g{fi')) in (2.9), and Cp^^nv = 

Cp. For cr"^, multiplication by 7^^ is just a Lorentz rotation, and does not 
change the anomalous dimension. Therefore, {q(T°'^'^^^Q^ , = (^cr"^75^^Q) 

and CaAli'.lj) = C^^^{ii\ii) [11], where a'f = |[72,7±]- 

There are 8 different Dirac matrices F (2.7) in 4-dimensional space. For our 
investigation of Cp, we can restrict ourselves to 

r = l, ^, 7^, 7°^, (2.14) 



because the other 4 matrices can be obtained from (2.14) by multiplying by 

The 1/m term in the expansion (2.8) was first investigated in [28], where the 
one-loop anomalous dimension matrix of dimension-4 operators Oj was found. 
The full one-loop corrections to for vector currents (and axial currents with 
anticommuting 75) were given in [29,30]. Some general properties of the match- 
ing coefficients Bi and the anomalous dimension matrix of Oj following from 
reparametrization invariance and equations of motion were established in [30], 
and the two- loop anomalous dimension matrix was calculated in [31,32]. 
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The sum in (2.8) includes the bilocal terms [33] originating from the insertion 
of the subleading terms of the Lagrangian (2.1) 

Cr I dxiT {j{0),Okix) + C^Omix)} . (2.15) 



The states \B>, \B*> in HQET matrix elements will be understood as the 
eigenstates of the leading HQET Hamiltonian with m = oo; they are m- 
independent. 

The terms in the sum with the derivative acting on the heavy-quark field can 
be obtained from reparametrization invariance [30]. Let T = 7^"^ • • -7""!. It 
can be decomposed into the parts commuting and anticommuting with ^: 

r = r+ + r_, r± = ^(r±^r^). 

The matrix element of the renormalized QCD current qVQ from the heavy- 
quark state with momentum mv to the light-quark state with momentum 
is 

I (Cr+ + Cr_) UgTu{mv) + i (Cr+ - Cr_) Ug^Tu{mv) . (2.16) 
In this equation, we may substitute v —>■ v + k/m: 



\ (Cr+ + Cr_) u^u{mv + k) 

+ \ (Cr+ - Cr_) [f + Vu{mv + k) 



Using u{mv -|- fc) = (1 + (2m)) u.u{k), we obtain the leading term (2.16) plus 
1 



4m 



(Cr+ + Cr_) UgF^u, + (Cr+ - Cr_) Ug m + 2|^r) 



Therefore, the qDhy terms in the sum in (2.8) are 

I (Cr+ + Cr_) qUpK + \ (Cr+ - Cr_) q {jVilp + 2%lpV) . (2.17) 



The coefficients of operators with the derivative acting on the hght-quark field 

are not determined by general considerations. These coefficients appear first 
at the one-loop level. We calculate the matrix element of the QCD current 
from the heavy quark with momentum mv to the light quark with momentum 
p (with — 0), expanded in p/m to the linear term, and equate it to the 
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corresponding HQET matrix element. In HQET, loop corrections contain no 
scale, and hence vanish (except, possibly, massive-quark loops, which first 
appear at the two-loop level). The QCD matrix element is proportional to 
u{p)T{p, mv)u{mv), where r(p, mv) is the bare proper vertex function. At one 
loop, it is given by Fig. 1. 



If we assume nothing about properties of F, then the term linear in p has the 
structure 



with LiX Ri = p- vlxl,p-v ^7^ x 7^, p ■ v •j^ju x 7*^7'^, 1 x ■^j^ x 7^1^. The 
coefficients Xi can be obtained, by solving a linear system, from the double 
traces of Dirac matrices to the left from F with Lj and those to the right from 
F with Rj, with Lj x Rj = x (1 + ^), -^pf x (1 + 1^)7/', 7^7^^ x (1 + ^)7''7'', 
ffx (1 + f)j)^ 7p|^ X (1 + f)'fri<^ . Now we can take these double traces of 
the integrand of Fig. 1, and express Xi via scalar integrals. Their numerators 
involve {k • p)"; putting k — {k • v)v + k± and averaging over k± directions in 
the {d — 1) -dimensional subspace, we can express them via the factors in the 
denominator. 

Now we assume 



mv k + mv \ k + p p 




Fig. 1. One-loop matching 



u{p) '^XiLiFRi u{mv) 



(2.18) 



^F = (7F|^ , (7 = ±1 , 7^F7'^ = 2ahr . 



(2.19) 



For our F matrices (2.7), 




(2.20) 



Then (2.18) becomes 



Xi + {x2 + 2x5) ■ 2h + Xs{2hy p ■ V u{p)Tu{mv) 
+ [x4 — X5 ■ 2h] u{p)T'^{mv) . 



(2.21) 



Performing the simple calculation, we obtain 
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u{p)r{p, mv)u{mv) 



2m ^(47r)<^/2- 



where 



u{p)ruy{o) 



(47r)<^/2 

b\p ■ vu{p)Tuy{0) + b2u{p)T;^Uy{0) 



(2.22) 



Cr = 



bl^-2 



bl 



[I - h){d - 2 + 2h) 

(d-2)(d-S) ' 
(d-2){d-8)-{d-5){d-A + 2h)h 

{d - 2){d - 2>){d - h) 
d-2-h 



{d-2){d-3) ' 



The zeroth-order coefficient cr has been found in [13,11]; the ffist-order coef- 
ficients for components of the vector current in [29,30]. 

For r = 1, the square bracket in (2.22) becomes b-p p-v u{p)ruy{0) , with bi = 
b\,b^ = 6^+262- We have checked these results by taking the trace of the whole 
integrand of Fig. 1 with (1+^)1^ and calculating the integrals. For F = 7°, 7"!^, 
the square bracket in (2.22) becomes 6r,iP • vu{p)'y'^luv{0) + br,2Pi_u{p)uy{0), 

with = br, ^^,2 = 2br, 6^,^,1 = bj^^ + 2bl^^, b,^^,2 = -261^''. We 
have checked these results by taking traces of the integrand with (1 + i))pai> 
and (1 + f)^ai>- An additional strong check is provided by the Ward identity: 
contracting the vertex function F"(p, mv) (for F = 7") with the momentum 
transfer {mv —p)a-i we obtain 

r"(p, mv){mv — p)a — mr(p, mv) + E(mi') , 



where F(p, mv) is the scalar vertex (for F = 1), and E is the heavy-quark 
self-energy. At the first order in p, this leads to 



61 - 6^ = 2 {c^^ - . 



(2.23) 



Our results (2.22) satisfy this requirement. 

In order to investigate the asymptotic behaviour of the perturbative series (2.12), 
we consider its Borel transform 



Sv{u) = 



u 



L-l 



(2.24) 



u=0 
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Formally one can invert this transformation and gets 

Cr = 1 + Jsriu) exp (- u] du . (2.25) 



However, if Sy{u) has singularities on the integration contour (which is the pos- 
itive u axis), then the integral (2.25) is not well-defined, and the series (2.12) 
is not Borel-summable. To deal with this singularities some prescription is 
needed, leading to ambiguities in Ct- 

In the large-/3o limit [11], 



Sr{u) = Cpl ^^'f:^^ f""^ bin - 2f - 2r){n - 2)u + 3u'' + u-5 

[ L [S U) 

_ 2(n-2)^-5 | 



(2.26) 



(the results for the components of the vector current were obtained in [7]). 
Expanding this Sr{u) in u reproduces leading large-/5o terms in (2.12) (in 
particular, in cf). This Borel image has IR renormalon poles at m > 0. The 
pole nearest to the origin (and thus giving the largest renormalon ambiguity) 
is situated at u — ^: 

St{u) — 1- (regular terms at u — ^) . (2.27) 

2 — u 



This leads to an ambiguity in the sum (2.25) of the series (2.12), the natural 
measure of which is the residue at the pole: 

(3o rn 



where is for rii flavours. This is commensurate with the 1/m corrections 
in (2.8). It is convenient to measure all such ambiguities in terms of the UV 
renormalon ambiguity of A [5] 

AA = -2C^e^/6^_ ^2.28) 



Then [11], 



ACr = -- 



2(n - 2)2 - 77(n - 2) - 



15" 



AA 



m 



(2.29) 
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In the following Sections we shall see that beyond the large-/9o limit S-r(u) has 
a branching point at n = | rather than the simple pole (2.27). We use the 
method based on renormalization group, which will be explained in Sect. 3. 



3 Spin-0 currents 



In the case of the currents with F = 1, ^, the leading term in the expan- 
sion (2.8) contains J = qh^, and the 1/m correction - four operators 



01 = qilpK = ida {qi°'K) , 

02 = q (^—iv ■ hy = —iv ■ d {qh^) , 

Oz^ij dxT{m,Ok{x)} , 
O^^ij dxT{j{Q),Om{x)} . 

Prom (2.15) and (2.17), 



(3.1) 



(3.2) 



The anomalous dimension matrix of the dimension-4 operators (3.1) has the 
structure 



7 + 



^ 


7^= 



(3.3) 



The operators Oi 2 are rcnormalizcd multiplicatively with 7 (this fixes the first 
two rows); the form of -81,3,4 (3.2) fixes the columns 1, 3, 4. The mixing of the 
local operator O2 into the bilocal operators 03,4 is described by [32] (these 
anomalous dimensions are called 7]^'° and 7^"^^ + 7"^^^ + 7"^^^ in this paper) 



a. 



V 9 3 / 



32 o 54 



160, 



, C'a TpTii 

9 J 9 



(3.4) 



The unknown coefficients i?2(A*')A*) F 
renormalization-group equations 



_dBl_ 

d log fjL 



1, ^ are obtained by solving the 



(3.5) 



with the initial conditions B^im^m) obtained by matching at /i' — /i — m. 
The ratio Bl{iJi' , jj) / Cy{h' , jj) does not depend on ji': 



B^{fi',fi) ]3l ''i^\''{as) das 



[ 



(3.6) 



(see (2.2)). The renormalization-group invariants B2 start from one loop: B 
^'2i«.(/^o)/(47r) + • • •. 



To find B^im, m) for F = 1, ^ with one- loop accuracy, we write down the sum 
in (2.8) via the bare operators: 



2 — 2f 

o Tn 

CF^^^r{e)br02o ± Oio + O30 + O40 

= ( Cpbr - ^£Mo2(m) + (other operators) . 

\ 2 y 47r£ 

Taking 7^ + 7™ = into account, we see that both br should vanish at £ = 0. 
The 0{e) terms of (2.22) give [29,30,34] 



Bl{m, m) = SCf^^ + • • • , Bt{m, m) = 12C^^^^ + 



47r 



47r 



and 



4^ = 8C^^ + 
47r 



, Bt^l2CF^^ + 



(3.7) 
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(using 7^ + 7^ = 0). 



An exact relation between and B2 can be derived. The QCD vector current 
and the scalar one are related by the equations of motion: 

idJS = id^r = rnojo = m{n')j{n') , (3.8) 



where 



is the n/-flavour MS running mass, and (2.11) 7^ = —GCpas/ (4:71) + •■■ is 
minus the mass anomalous dimension. We separate = (j ■ f )f " + and 
substitute the expansions (2.8) with (3.1), (3.2). The matrix element of (3.8) 
from the heavy quark with momentum mv to the on-shell light quark with 
momentum p reads 



2m 



m(//')Ci(//',/x) jl + — 



p ■ V + r 



1 \ p ■ V + r 



where 



<q\iJdxT{j{tJ,), {Ok + CMOM)J \Q>-Pa <q\r{lj)\Q> 

<q\j{i^)\Q> 



T — Ti'^^v At the leading order in 1/m, this yields [11] 

m _Ci(/i',/i) ^^C^ ^gg^ 



At the first order, we obtain 



Ci{iJ,',fj,) C^{n',n) \C^{fj.',fj,) 



El 



Bl 



or ^ - ^ - 2 
Ci C4 



, C4 



- 1 



(3.10) 
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Note that (2.23) is just the one-loop case of this general result. The one-loop 
results (3.7), of course, satisfy this requirement. 

As discussed above, these results don't change if we replace q — > 57^*-'. Now 
the leading term is j" = Ql^^h^, and the definitions of Oj (3.1) are changed 
accordingly. We define the matrix elements 

<0| {qit^'Q)^, \B> - -irriBfM) , 

<mit''rQ\B> = -tfBP% , (3.11) 
where 



/I 



/I 



.<(/^o) 




(3.12) 



The HQET matrix elements are [35] 



<0\j{fi)\B> ^ -i^Fin) , 

<0\Oi(iJ,)\B> = -<0\O2\B> = -i^/^AF{n) , 

<0\Oz{ix)\B> = -i^/^F{fi)GM , 



(3.13) 



(in [35], Gk and Gm are called 2Gi and I2G2), where A = — m is the 
S-meson residual energy, 

F{^) = F I^Y''' K,{a,{n)) , (3.14) 



and F is /^-independent and is thus just a (non-perturbative) number times 
A^. The hadronic parameters Gk,m{l^) obey the renormalization-group equa- 
tions 



— = 7 )A, 

alog^u 

^§^ + lra{as{^))GM = r{a,{^))K. (3.15) 
a log// 

Their solution is 
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2/3o ^ 47r ^ 



as (ft) 



^■(q,.) 7^ \ das 



2(3{as) 2PoJ as 



All 



X 



- A 



7™ (a 



as\ 2/3o da^ 



47r 



a.. 



(3.16) 



where Gk and are again /^-independent and thus are just some (non- 
perturbative) numbers times A^. 



Taking the matrix element of (2.8), we obtain 



fB 



C^(/^^/x)F(/x) 



X 



1 + IT- {cl{^J)X + GM + Cm{^J)Gm{^J)) 



2m 

where F = 1, ^, and 



(3.17) 



, Ci{t,) = l-2 



Substituting the solutions of the renormalization-group equations, we arrive 
at the explicitly //-independent expressions 



fB 



X 



'as{fio)\^ CrF 



l + —\C'j,A + Gk + CmGr, 



47r ] 



where 



(3.18) 



Cl^\-2\ 







Bl 



Cr 

2/3o 47r J [2P{as) 



7o ^ 



) 2/3o; as 



15 



J 2p[as) as 



At the next-to-leading order, 



fE 

fB 



+ 



2m 



2A 







47r 
47r j 



7mO 



47r 



47r 



A 



(3.19) 



where 



= (1 T 1) (cf - cr) 

2/?o (7mo - 2/3o) 
(here 7o^ + 7^ = 0). 

In the large /3o limit, the leading short-distance coefficient C*r in (3.19) has the 
IR renormalon ambiguity (2.29). If we change the prescription for calculating 
the integral (2.25) near u — |, we should adjust the hadronic parameters A, 
Gki Gm accordingly. The UV renormalon ambiguity of A is given by (2.28) in 
this hmit. The UV renormalon ambiguities of Gk,m{l^) [7] 

AGkifi) = ~AA , AGmifi) = 2AA (3.20) 
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7mO 



7r + 7r 

2/3o 



+ 



Pi (to' + 1^) 



are /^-independent. We can obtain them by a direct calculation, performing a 
similar analysis as in [10], see Appendix A. In the large /So limit, we obtain 
from (3.16) and (3.20) 

a4 = --AA, AG'^= 2--^ AA. (3.21) 

2 V 7mO / 

In the 1/m expansion for the meson decay constants (3.19), F = 1 and f 
correspond to n = 0, ?7 = —1 and n = 1, = 1 in (2.29), and the IR renor- 
malon ambiguities (2.29) are — (3/4)(AA/m) and (l/4)(AA/m). They cancel 
with the UV renormalon ambiguities of the hadronic matrix elements (2.28), 
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(3.21) in the 1/m correction in (3.19). In fact, the results (3.20) were first 
obtained [7] from the requirement of such cancellation for F = 7_l, by solv- 
ing a system of two linear equations, and later confirmed [11] by considering 
all possible V (this gives three equations, thus providing a consistency check). 
Here we demonstrate the cancellation of renormalon ambiguities by a direct 
calculation. 

This cancellation should hold also beyond the large /3o limit. By dimensional 
arguments, the UV renormalon ambiguities of the //-independent hadronic 
parameters A, G^, Gm must be equal to h.y^s times some numbers: 



AA = 7VoAo, AG'fe = --AriAo, AG^ = 2 - ^ Aq , 

2 V 7mO / 

Ao = -2C^e^/6^ . (3.22) 

The normalization factors are unity in the large /?o limit, — 1 + C(l//5o); 
in general, they are just some unknown numbers of order one. Using 



^MS = /^o exp 
K{as) — exp 



27r 



/5oas(/Wo) 
^ 1 



V 47r ) 

27r ^ A 



_9 

2/3; 



° K{as{iJ.o)) , 



(3.23) 



^2/3(a,) /3oa. 2131) 



, _ /3o/32 - 131 as 

47r + 



/xo = e-^/^m (2.5), we can represent the UV renormalon ambiguity of the 
1/m correction in (3.19) as exp[— 27r/(/3oQ;5(//o))] times a sum of terms with 
different fractional powers of Q;s(//o)/(47r). ^ In order to cancel this ambiguity, 
we should have the branching point 



Sr{u) = E 



(i 



(3.24) 



instead of a simple pole (2.27) (here S^^{u) is regular ect u = We define 
the IR renormalon ambiguity of Cr - generalizing the prescription to take the 
residue of the pole - to be the integral of (3.24) around the cut divided by 
2ni: 

^ It is convenient to replace log[Q;5(//o)/(47r)] — > [{ag^/J^o) / i'^'^))'' ~ l]/<^) a-^d take 
the limit 5 ^ at the end of calculation. 



ACr 



exp 



27r 
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The requirement of cancellation of the ambiguities gives 



Sr{u) 



1+ 



01 

2^ 



i- / 1 

I log 



2/3o 
3 



^P{1 + 



7mO 



7mO 



(3.26) 



where 



/3l 



r 1 + 



2/32 2/?o. 



/5l TmO 



-'ml 



2^ 



/?o(2c 

ml 

/3l - PolmO 



(3.27) 



In the large Po limit, the formula (3.26) reproduces the known results (2.27), 
(2.29). 

The asymptotics of c£ (2.12) at L » 1 is determined by the renormalon 
singularity closest to the origin (see (2.24)). At n ^ 1, 



T./ -.N , a f. a{a + l) a(a^ - l){3a + 2) 
r(n + a + 1) = n! 1 + + — ^f^^ + 



2n 



24n2 



and we arrive at 



cl^, = 2Cf n\ (2/3o)" (2/3on)^ + + . . .j 



7o 



7o™ 



+ 



J// 



log2/?on±l+ , 
2po 7mO 2/^Jon 



+ 



2 V 7mO, 



N2{2pon) W 1 + 



'-"ml 

2/3ori 



+ 



(3.28) 



where 



r 2/?o/?2 - 2/5o'A - 3/5? _ r 



, c^i — + 7q 



2/3o^ ' 
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n _ _ 2(/3i+/3g)-/3o7mO 
C-ml ~ '^rnl 7mO , n2 



(3.29) 



The result (3.28) is model-independent and the powers of n are exact. However, 
the normalization factors Ni cannot be determined within this approach. The 
leading term at n ^ oo formally is the logarithmic term, because 7^0 is 
positive (sec (2.4)). At moderate values of n, all leading terms are of similar 
importance. The function Sy'^{u) in (3.24) has singularities at m = 1 (IR) and 
U——1 (UV), and thus gives exponentially smaller contributions with (i/^o)" 
instead of (2/3o)". 

The matrix element of (3.8) is 



vriB 



m 



(3.30) 



Substituting (3.17) (or (3.18)) and using the relation (3.10), we obtain 

/!(//) Ci(/x',/x) 



B 



1 + 

m , 



(3.31) 



The ratio of the quark masses is given by (3.9). Naturally, it contains no l/m 
corrections with i?-meson hadronic parameters; it is just a series in as{no), 
see the Appendix B. The Borel image of this series is 



1+ 



I3l 
2^ 



- f3l 



Pi 



(m/m)/ 
C2 



2/?o' J ' 

{m/rh) 



(3.32) 



2/3o^/33 - m + A (A - 2/32)(2/3o/32 - 0{\ 



and the asymptotics of c^^"*^ at L ^ 1 is 



(m/m)// 
'^2 



2Pon {2Pony 



+ 



{m/m)» _ (m/m) 2/3q/?2 — 2/3q/5i — SfS^ 
Ci — Ci 



(3.33) 
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{m/m)>f {m/rh) {m/m) 

24(3^(3s - 12131131 + - 2/32)(36/?o/?2 - 27/?^ - IQ^lP^ - sp^) 



These results are equivalent to [9] . 

We can try to estimate the unknown normalization constant A'o following [36,37] . 
The function 

1+ At ~ pW™) / „ , \ 



still has a singularity at m = ^ due to S^^^{u) in (3.24), but has a finite limit 
at — > I — 0. The radius of convergence of its expansion in u is thus |, but 

the series should converge aX u — |. Therefore, we can calculate Sm/m (|), 

and hence Nq (3.32), from this expansion. Substituting c^'"/™) for L < 3 from 
Appendix B, we find 

Nq = 0.288 ■ (1 + 0.075 + 0.630 + •••)- 0.491 . (3.34) 



The three-loop correction turns out to be large, which casts some doubt on 
this estimate of A^o- 



4 Spin-1 currents 



In the case of the currents with F = 7", 7"^, the leading term in the expan- 
sion (2.8) contains J" = q'j^^hy, and the 1/m correction - six operators 



Ot = q{-iDl) K, 

02 = q (-IV ■ 72 /i„ = -tv-d {q-ilK) , 
Ot^ij dxT{r{Q),Ok{x)) , 

O^^ij dxT{r{Q),0^{x)} . (4.1) 
Note that 

0« - Of = tdl{qK) . (4.2) 
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Prom (2.15) and (2.17), 



BJ^ = 2C^ , -B]^ 



Br 



a 



7± ' 



B^ — CmC'yj_ 1 



-2-Dl 
-^6 



B 



B\ 



a 



■y±f ' 



(4.3) 



The anomalous dimension matrix of the dimension-4 operators (4.1) has the 
structure 



7 + 



(OO -fa 76 0\ 



7, 7b 

7^ 
VO 7^72^^0 7^; 



(4.4) 



The operators 02,4 renormahzed multiphcatively with 7, which determines 

the second and the fourth row. The same holds for Of — Of (4.2), and hence the 
first and the third rows coincide. Furthermore, the form of -81,2,5,6 (see (4.3)) 
fixes the columns 1, 2, 5, 6. Renormalization of the operators (4.1) is discussed 
in Appendix C in detail. 

The unknown coefficients B^^{fi',fi) for F = 7^, 7^^ are obtained by solving 
the renormalization-group equations 



dBl 



= (7 + 7a) Bl + 7asr + 7r^6 , 



dlogfx 

- 7^1 + 76 (Sf + Bl) + ^^Bl + 7-S, 



dlogii 



(4.5) 



(where (2.11) '^'^^.^ = ICpagj (47r) + ■ ■ ■) with the initial conditions B\^{m^ m) 
obtained by matching. The ratios Bf (/x', /i)/C'r (//■', yu) don't depend on /x': 



Bl{fJ^',fi) _ ( as{n) 



TgQ 
2/3o 



C'r(/i',/i) \Ois{lJio)) 



K 



laMlJ')) 



Bl 



as(Mo) 



P[Ois) dig 



(4.6) 
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as(Mo) 



2/3(«,) 



Cr(/x',/x) Cr / 2/3(a,) 



+ 2^ 



-1 



do;. 



TiGts) das 
2f5{as) as 



o:s(/io) 

- C^as{iio)'^ 



2/?(q;,) 



2/3o 



(4.7) 



(in the last formula, the running {/i' , /i) / Cr{lJ'' , IJ>) corresponding to the 
integration variable as is understood). 

To find ^{ni, rn) for V = j±, j±'^ with the one-loop accuracy, we write down 
the sum in (2.8) via the bare operators: 



(47r)'^/2 

A-KS 



+ {cFhv,! T 760 - ^° OJ(m) + (other operators) 



The values of 6r,i at e = have to cancel these anomalous dimensions. The 
0{e) terms of (2.22) give 



Bl^ (m, m) = 4Cf + • • • , (m, m) = -4(7^.- 



47r 

Q:s(m) 



47r 



+ 



53^^''(m, m) = 2Cf^^ + • • • , 52"''(m, m) = -QCf^^ + 



47r 



47r 



The results for F = 7j^ were obtained in [29,30]; those for F = are new. 
Using the one-loop anomalous dimensions, we get 



3 — 



2^ Osipo). 



26 a^(/^o) , 



57^ 



«^ «£(/f0) , 



A7X)< 2 oJ^IIq) 



(4.8) 



We define the matrix elements 
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<0\q^''Q\B*> = mB./B*e" , 

<0| {qa'^^Q)^, \B*> = (e>|. - e^p^.) , (4.9) 

where is the S*-meson polarization vector. 
The HQET matrix elements are [35] 

<0|r(/i)|fi*> = v/^i^(/i)e", (4.10) 
<0\O^{^)\B*> = <0|O^(/i)|5*> = -7^AF(/i)e". 

The matrix elements of and O3 are equal, due to (4.2). However, the 
formulae [35] for these matrix elements hold only at the leading order. Let's 
define 

<0|O?(/x)|S*> = <0\O^(ii)\B*> = ~^/?^AF(ii)R(as(i^))e^ , (4.11) 

where R = 1 + 0{a^). It obeys the renormalization-group equation 

dR 



dlog/i 



+ 7„i? + 376 = 0. (4.12) 



Following [35], we define 

<0\Ol{^i)\M> ^^^Tr^lFM, (4.13) 



where Of is defined in (C.l), M = B or B* , and is the corresponding Dirac 
structure. If we take F = 7aF', then 

<0\Ol{ii)\M> = ^F2(//) TrF'M . 



Taking into account (C.8) and 

<0|Oi(A()|M> = <0|O^(/i)|M> = -^AF(//)TrF'M, 



we obtain at the next-to-leading order 

1 - la 
FM = --AF(/x)i?(a,(//)) , R{a,) = 1 + -7ao^ + C(«D ■ (4-14) 
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The general result for R{as) can be derived by solving (4.12) and requiring the 
absence of fractional powers of ag (or by requiring that (4.14) is reproduced): 



if 7/37^ 



a. 



.(4.15) 



The matrix element of is [35] 

<0\Ot{l^)\B*> = V^^(/^)G'ik(/^)e". 



(4.16) 



However, the formulae [35] for the matrix element of Og hold only at the 
leading order. Following [35], we define 



<0|O3™(/x)|M> = ^F(/.)G(/x) Trr^^a,,Ma^^ , 



(4.17) 



where is defined in (C.9). The 5-meson matrix element G^(/x) (3.13) is, 
at the next-to- leading order (C.14), 



(4.18) 



For B* we define 



<0\O^\B*> = -- [GM + Rmiasii^))A\ V^^(/^)e" , 



(4.19) 



where at the next-to- leading order (C.18) 



(4.20) 



It obeys the renormalization-group equation 

+ ^mRm + 7™ + iTR + 372™ = . 



dlogfj, 



(4.21) 



Its solution (which contains no fractional powers of and reproduces (4.20)) 
is 



7mO 
, 2/3o 



^'^^{as)+^TMRM + 3^^{as) 



(4.22) 



X 



/ 



2/3(a,) 



a. 
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Taking the matrix element of (2.8), we obtain 



X 



2m 

where F = 7^, 7^^, and 



(4.23) 



1 



Cr(/i', /i) 



Substituting the solutions of the renormalization-group equations, we arrive 
at the //-independent expressions 



B* 



fT 
J B" 



TO 

2/3o CrF 



X 



1 + — ( C^-^ + Gk — ^CmGrn ^ 



7mO ■ 

2/3o 



2m 



47r 



where 



(4.24) 



CAM 

3 -1 



> + 1 i?(a.(/.)) 



Cr(//',//) 



2Po 47r y y2(3{a,) 2Po J 3 



7''("s) 7o \ das 



(4.25) 
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At first sight, it is not obvious that does not depend on fj,. However, 
differentiating it in log/i and taking into account the renormahzation-group 
equations (4.5), (4.12), (4.21), we obtain zero. Therefore, the expression with 
II = can be used: 



Bl_l 
Cr 3 



M/fo) 

2/3o ^ 47r 



I C^/ C, 

1 - 

f / tM 7o 



+ 7r^-7a(«s(/^o)) 

1 I C'r 

'y'^iois) 7o \ das 



a. 



-CmK_^^{as{i^o))Rm{as{l^o)) 



Os(Mo) 

+ |Cm«,(/Xo)W J 




2m0 



2p(Q;s) 



(4.26) 



Comparing this with (4.25), we can rewrite B^{ii', n) / Cr(/^', yw) in a form which 
seems different from (4.7) but is equal to it: 



S[(^^ S[ _1 
Cr(/i',/i) C'r 3 



7X , ^ ^Mi^ , ^^^^^^^^ 



Cr(/i',/i) 



/ 



+ 



-1 



B. 



3_ 



» + —K_^^{as{i^o)) I i?(a3(/xo)) 



'7^("s) 7o \ 



V 

-^log , . 

1 1 

- gC'^(/^)-Rm(Q;s(/^)) + :^CrriK^^^{as{iio))Rm{as{fJ'o)) 



as (mo) 



(4.27) 



(to convince oneself that this is equivalent to (4.7), one can check that they 
coincide at /j, — /iq, and that (4.27) obeys the renormalization-group equa- 
tion (4.5).) 



At the next-to-leading order. 
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a,(/io)\^ CrF 



fl 



An 



+ 



2m [\ 
+ Gk 



Wo ^ + 3 I " ^ J + '^^^^ + 



1 



1 + Cml — : V 



47r 



47r 



A 



(4.28) 



where 



= 1(1 ± 1) (cr - c?) - \hl, - h\ - i^c^i T i7ao + + |72™0 

^ ^ 7mO 

2/3o 2/32 6/3o (7mo - 2/3o) ' 

In the largc-/?o hmit, the IR renormalon ambiguities ACr (2.29) for F = 
7±; 7±^ (having n = 1, 77 = — 1 and n — 2, rj — 1) are (ll/12)(AA/m) 
and (5/4)(AA/m). They cancel with the UV renormalon ambiguities of the 
hadronic matrix elements (2.28), (3.21) in the 1/m correction in (4.28). 

The Borel images are 



Sr{u) 



1+ 



2^ 



3 V 7mO, 



^7Vi-l 2-^ 
2 ' 3 V 7mO, 



l + c[' {\-u) + 



2/9o^ 



+ Cai(|-«) + 



^0 



TmO 
2/3o 



+ 



,(4.29) 



(see (3.27)), and the perturbative coefficients at n ^ 1 are 



cl^,^2CFn\ (2/3o)"(2/3on)^ 1 + 



2i%n 



^ -J V ImO, 



7o 



7„o / 2/3on 



1 + 



'-ml 



(4.30) 
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(see (3.29)). 

The ratio /b*/ fs* at the leading order in 1/m is given by the perturbative 
series in Q;s(/io) with 



(/r* / Ib* ) /-I 



, (4.31) 



Cf\ - + ^)fi'o + (2C3 + |7rMog2 - f TT^ + ^) Cp 



+ (-4(3 - 1-^ log 2 + f - 1) + (-1^^ + f ) T, 
+ [(32C3 - f ) + (-78C3 + ^) C^C^ + (42C3 - #) C; 

250-<-< rp I 80/^ 1 1 
-I Po 

I /'625/-<3 I 461 /-(2 /-I , 2375-^ -^2 511-^3 

I 220-^ m I 140,^2 \ 1 
+ -g-<^F<^A-LF + —<^A-^f)-^ 



Ca 



- I (75C| + 25CfCa + 21C1) {UCp + 7C^) 
+ iC^(llC^ + 7C^)^^| 



(from the result in [11], omitting the rric ^ effect, and the three-loop anoma- 
lous dimension 7^ of the tensor current [24]). This ratio is, from (4.24), 



Ib* 



ifg* I fa* ) 

Cm 



l-Afi+eif/^-)^ + 
3m \ 47r 

Cm , 



Therefore, the Borel image of the perturbative series is 



(4.32) 



2 Cf% 



3/1 



1+ 



/3l 
2^ 



l + cS^-/^-^'(i-«) + 



Cl 



Pi 



Cai + Cl 



(4.33) 



and the asymptotics of Q^B*/fB*) ^ ^ ]^ jg 



28 



/5l 



n+1 = —CfN^uI (2/?o)" (2/?on)^ | 1 + ^ 



Cl 



2/3on 



+ 



, (4.34) 



5 Results and conclusion 



Our main result is for the ratio of two (in principle) measurable quantities, 
/b*//b- It is given by the perturbative series in as{iio) with 



„{fB*/fB) _ 

C2 — 



(5.1) 



-3/3^ + (-8C3 + f vr^ log 2 - f TT^ + f ) 
+ (4C3 - log 2 + fTT^ - 6) + (f TT^ - f^) T, " 

(from [11], omitting the mc 7^ effect). This ratio is, from (3.19) and (4.28), 



a 



7x 



^ 2 
3m 



7mO 



47r 



Air 



where 



A 



(5.2) 



JfB*/fB) _ ^ (^1± _J \ 

^6Cf 

The Borel image of the perturbative series is 



_ 4Cf 8(47r^ + 3)Cf - (Stt^ - 93)Ca 
~ 3 ^ 27(/3o - Ca) , 



SfB*/fB{u) = 2 



2/3^ L 



1 



- 2- 



7mO, 



a-) 

iv^i-«)^(i + cL^r/^-)'(i -«) + ...) 



(5.3) 



where 
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I3l - PolmO 



Mb*Ub)i _ MB*/fB) _ - Pi 
The asymptotics of the coefficients is 



MB'Ifs) 
-71+ 1 



- 2- 



-Cpnl (2/?o)"(2/3on)^ 



MB*/fB)" 

1_^ + ^A1 + 

7^0 2/3on 



(/s*//s)" 



^)7V,(2/3on)-^ 1 + + 

7mo y \ 2/?on 



(5.4) 



where 



Cm 



7mO, 



(/s*//s)" _ ^ , J/s*//s)' , (A - /3o7mo)(/3l + 2/3o - /?o7mo) 



-'ml 



— CjTil + Ci 



+ 



Substituting the parameters, we have the perturbative series 



fs 3 TT 



--C3 + ;^7rMog2 + -7r2 + — ) 



1. 2 
-Cs + — • 
9^ 27 



36 ; V TT 



(5.5) 



The asymptotics (5.4) becomes 



81 V 3 y V 3 / \ [ 25 V 250000/ n J " 
2,. /50 / 40157 1 \1 

Definite quantitative predictions cannot be made, because the normahzation 
coefficients Nq^2 are unknown. To have some idea about the growth of the 
coefficients, we present them in Table 1. The coefficients cl/4^ oias{iiQ)/'n are 
given in three columns. The first column shows the exactly known ones (5.5). 



30 



The second column shows the results of the large-/3o limit. The Borel image 
Sfg^/fg{u) in this limit is, from (2.26), 

_ r(i + ^)r(i-2^) 
Sf../fM - -^c, 

Expanding it at m = (2.24), we get 

3C^(8C3 + 27r2 + f)/5o', 

(144C3 + f + 287r2 + 93) /?o\ 

Cf (144OC5 + 16071^3 + 840(3 + 367r^ + ISOtt' + ^) /?o' , ■ ■ ■ 

This limit reproduces c^^'^*^^^^ and the /3o-term of 02'^*^^'^'' (5.1). Finally, the 
third column shows the asymptotics (5.6) at Nq = N2 = 1. Let's stress once 
more that this is not the result of QCD, but simply a numerical illustration 
of the typical behaviour of the perturbative coefficients. The large-/9o result 
includes not just the pole at -u = |, but the whole function 5'/^,//g('u) (5.7). In 
contrast to this, the asymptotics (5.6) is determined by the nearest singularity 
&X, u— ^ only, but includes all powers of /3o, not just the highest one. To show 
the rate of convergence of the 1/n expansion (5.6), the asymptotic result is 
also expressed via the leading term and the 1/n correction. One can see that 
the accuracy of our next-to-leading order results at L < 10 is not high. Finally, 
the last column shows the complete L-loop contribution to according 
to (5.6) with ATq = iVa = 1. We use the value as{e-^'^mb) = 0.299 obtained 
by RunDec [38] . The smallest contribution seems to be the 3- loop one, and it 
is about 4% (though this small value is due to a partial cancellation between 
the leading order and the next-to-leading one, and thus is not quite reliable). 
Therefore, calculation of this 3-loop correction is meaningful (and it is actually 
possible, using the technique of [39,40]), while there would be no sense in the 
4-loop calculation. 

For the ratio m/rh, we obtain (from the Appendix B) 



(5.7) 



lfB*/fB) 
C4 

MB*/fB) 
MB*/fB) 



— 1 _L 891 Q'g(/^o) I ( 173 I 1^2 i^„r, I 1^2 I 168550145^ / ^sjl^'o) 

- + 1058 ^ + l~138^3 + gTT lOg ^ + gTT + 40297104 ) I ^ 

_L 1^ 188^ I 50225 A 1439 „2 A 4396763 /- 47 i _„4 r, 14 2 i 2 r, 
+ (,-^«4 + l96^C5 - l32-7r Cs - 138012-C3 " 162 ^ - gjTT lOg 2 

_ 402485 2 1 r, I 461 4 , 220317449 2 , 17341442069927 A / ^s{^o) \ 
85698 " ^ 7776 " 20567520 " 767418048576 j 



TT 
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Table 1 



The perturbative series for fs* / fs 





_JjB*/fB) 


\jKjl.Lyjm 


exact 


large (3q 


asymptotic 


1 


0.67 


0.67 








2 


4.06 


2.08 








3 




29.17 


47.26 = 95 • (1 


-0.501) 


0.0406 


4 




333.26 


902.44 = 1363 • (1 


- 0.338) 


0.0736 


5 




6342.19 


18699.42 = 25103 • (1 


- 0.255) 


0.1450 


6 




128998.30 


449431.53 = 565332 • (1 


- 0.205) 


0.3311 



+ ... + o(^\. (5.8) 
\ rrih J 

The asymptotics (3.33) becomes 



(m/m) 16 . 

c„+i = y ^0 n\ 



+ 



880261 
7187500 



50 

y 

C3 + 



50 



-n 



231 
625 



1 + 



688161953 1 
1653125000 n 



(5.9) 



4 2,^ 4 2 8332134087653830381 \ 1 

TT^ log 2 H TT^ H — 

625 ^ 625 491908007812500000007 v? 



It depends on just one normalization constant A^o- In the large-/3o limit, from (2.26) 
we obtain [5] 



m/ih 



U) 



6Cp 



T(u)T(l 



1 



r(3-u) ^ 2u 



(5.10) 



and 



4"^/-) = (72OCI + IO8OC5 + 1207r^C3 + 9OC3 + + 97r^ + fn' + f) . . . 

(the terms with the highest powers of in the Appendix B are reproduced). 
Numerical results are shown in Table 2. For L = 3, the 1/(L — 1) expan- 
sion (5.9) seems to converge well; comparison with the exact 3-loop result 
from Appendix B suggests that the normalization factor A^o is smaller that its 
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large-/3o value 1, namely, Nq ^ 0.27. This conclusion is in a qualitative agree- 
ment with the estimate (3.34), especially if we omit the problematic 3-loop 
correction in it. 

Table 2 



The perturbative series for m/rfi 



L 




cent. 


exact 


large Po 


asymptotic 


1 


0.84 


0.5 






2 


4.53 


7.37 






3 


56.37 


36.23 


209.88 = 169 • (1 + 0.208 + 0.032) 


0.180 


4 




641.71 


2833.12 = 2457- (1 + 0.139 + 0.014) 


0.231 


5 




9062.48 


50650.18 = 45543- (1 + 0.104 + 0.008) 


0.393 


6 




206941.30 


1121489.20 = 1030382 • (1 + 0.083 + 0.005) 


0.826 



Finally, our last ratio (4.31) is 



fs* _ 


-, 707 "^(/^o) 




3174 ^ 


+ • 






\ rrib J 



-(^^C3 + ^7rMog2 + A,^ 



102436609 \ 
120891312 J 



(5.11) 



The asymptotics (4.34) becomes 

c«-/'-) = -«7V„«! (50y/50 xiS/ 792338047 1 ^ 

9°V3;V3;V 165312500071 J ^ ' 



It also depends on just one normalization constant A^o- In the large-/5o limit, 
from (2.26), S jt /fg,{'u) differs from (5.10) by the factor — | (this reproduces 
the terms with the highest powers of Po in (4.34)). Numerical results are shown 
in Table 3. 

If we neglect subleading as corrections, then, from (4.32) and (3.31), 

fl/fB^ = {fE/fB)''^' ■ (5.13) 



This equality also holds at the first order in l/(3o, to all orders of a^. There- 
fore, the ratio of the perturbative coefficients (5.12) and (5.9) is — |, up to 
corrections suppressed by 1/n and 1/Pq. Similarly, if we neglect subleading 
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Table 3 

The perturbative series for f^*/ fB* 



Lj 




cont . 


exact 


large /3o 


asymptotic 


1 


0.22 


0.17 






2 


1.04 


2.46 






3 




12.08 


42.88 = 56 • (1 - 0.240) 


0.0368 


4 




213.90 


688.14 = 819 -(1 -0.160) 


0.0561 


5 




3020.83 


13361.95 = 15181 • (1 - 0.120) 


0.1036 


6 




68980.43 


310536.81 = 343461 • (1 - 0.096) 


0.2288 



corrections, including those suppressed by [Q;s/(47r)]'''"'°/^^^°\ then, from (5.2) 
and (3.31), 

Therefore, the leading asymptotics of the perturbative series for fs*/ fs (5.6) 
and m/m (5.9) are related by 

/bV/b = (m/rn)-" . (5.14) 

The term with N2 in (5.6) violating this relation is suppressed not only by 
(2/?oT^)^^^^^, but also by a small numerical factor |. This approximate relation 
was first noted empirically at the 2-loop level in [11], with the exponent a = |, 
which is very close to 14/27. 

Let's summarize our main results. 

(1) The behaviour of the Borel images of perturbative series near the leading 
singularity u — ^ior the matching coefficients (3.26), (4.29), and for the 
ratios m/m (3.32), f^*/ fs* (4.33), and fs*/ fs (5.3) has been found. The 

powers of i — M are exact; further corrections are suppressed by positive 
integer powers of | — -u. The normalization factors A'"o,i,2 cannot be found 
within this approach; they are some unknown numbers of order unity. 
Logarithmic branching is a new feature of this problem; it follows from 
the fact that the anomalous dimensions matrices cannot be diagonalized. 

(2) Asymptotics of perturbative coefficients at L ^ 1 for the matching 
coefficients (3.28), (4.30), and for the same ratios (3.33), (4.34), (5.4) have 
been found. The powers of n = L — 1 are exact; further corrections are 
suppressed by positive integer powers of 1/n. Logarithmic terms follow 
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from the same property of the anomalous dimensions. 

(3) The coefficients and B2 of the subleading operator O2 for the currents 
with r = 1 and ^ are related by (3.10). One-loop results for the generic 
r (2.22) and for the tensor current (4.8) are also new. 

(4) The heavy-quark symmetry relations [35] for B- and S*-meson matrix ele- 
ments of subleading operators get non-trivial radiative corrections (4.14), 
(4.20). 

Note added. The three-loop anomalous dimension of the heavy-light quark 
current 72 has been calculated recently [48]. Therefore, the coefficients C2 
in (2.12) are now known, for all F. 
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A UV renormalon ambiguities in the large-/3o limit 

Ultraviolet contributions to the matrix elements of 03,4 are independent of the 
external states, and we may use quark instead of hadron states (see (3.13)). By 
dimensional analysis, the UV renormalon ambiguities of the matrix elements of 
O3 4 are proportional to AA times the matrix element of the lower-dimensional 
operator j with the same external states. We consider transition from an off- 
shell heavy quark with residual energy < to a light quark with zero 
momentum, this is enough to ensures the absence of IR divergences. For O3, 
all loop corrections to the vertex function (see Fig. A.l) vanish. The kinetic- 
energy vertices contain no Dirac matrices, and we may take | of the trace on 
the light-quark line; this yields at the vertex, and the gluon propagator 
with insertions is transverse. There is one more contribution [10], which we 
have to take into account. The matrix element -F of J should be multiplied by 

1 /2 

the heavy-quark wave-function renormalization , which contains a kinetic- 
energy contribution. This contribution is known to have an UV renormalon 
ambiguity [10] 

3AA 
2 m 

thus giving — (3/4)(AA/m)F as the ambiguity of the matrix element of O3. 
This must be equal to F • AGk/ (2m), and we recover the first result in (3.20). 
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Fig. A.l. Matrix element of O3; renormalon chains are inserted into the gluon prop- 
agators 

For O4, a straightforward calculation of the diagram similar to the first one in 
Fig. A.l gives the bare matrix element of the usual form (see [11,10]) 



1 - Fje, Le) 



with j3 = /5oas(/^)/(47r), 



F{s,u) = -2drCFC^^- 

m 



-2uj) 



1u 



.uV{-\^2u)T{\-u) 



V{2 + u-£) 



where 



D{e) = 6e^"r(l + €)B{2 -e,2-s) 



l + -e + 



dr = - 



1 + + 



1 + 



2 2 



;i - 2r;(n - 2))^ - 3 



The renormalization-group invariant matrix element has the form (2.25) with 
Ho = -2uje~^l^ and 



S{u) 



F(0,ii)-F(0,0) 



u 



-2drCFCm{-2u) 



LU 



m 



fl=HO 

r(-i + 2M)r(i-M) 
r(2 + u) 



+ 



2u 



Taking the residue at the pole -u = |, we find the UV renormalon ambiguity 
{dr / 3)Cm{—2Lj) {AA/ m) times the matrix element of J. For 5-meson, dr = 3; 
comparing with (3.17) gives the second result in (3.20). 
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B The perturbative series for m/rh 



The ratio of the on-shell mass m and the renormahzation-group invariant MS 
mass m is the series 



m 



1 I \ " „im/m) 



oj 



in the nrflavour Q;s(/io) (a*o = e~^/^m). Using the relation [41,39]'^ between 
m and fh{m) (omitting the rric effect known at two loops [41]) together with 
the 4- loop /3-function [43] and the mass anomalous dimension [44,45] , as well 
as [46,47] 



l + ^(32C^-39C^)r^[^ 



+ 



we obtain 



(m/m) 



(m/m) 
'-•2 



I - (f + 8Ca) ^ + 3 (IIC^ + 7Ca) §i 
Po Po 

<^f| (tt' + f)(3'o + (-6C3 - Stt^ log 2 + Stt^ + 



+ (12C3 + 47r^ log 2-5n' + ^) + (47r2 - f ) 



+ 



+ (66C3 - '-f) CfCa + (-66C3 + 1^) 
, 1 



+ -^(l350C| + 504C2^A - 3948Ci.C2 - 483C; 



3520CfCaTf - 2240ClrF 



'/2 



+ I (-15C| - 5C^C^ + ICl) {llCp + 7Ca) ^,3 

Po 

+ |C^(11C^ + 7C^)^^|, 



+ 



(l28a4 + IOOC3 + f log^ 2 + f tt' log' 2 - 327r' log 2 



^ The three-loop coefficient in it had been found numerically [42] before the ana- 
lytical result [39] was obtained. 
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22^4 I 47^2 I 383\r^ 

+ (-64a4 - 6IC3 - I log^ 2 - f TT^ log^ 2 + IGtt^ log 2 

+ (48C3-K-f)T^]/?^ 

+ (768a4 + 8OC5 + 47r^C3 + 305(3 + 32 log^ 2 - 327r^ log^ 2 

- 5087r2 log 2 + Itt^ + ^TT^ + ^) C| 
+ (-384a4 - 2OOC5 + 767r=^C3 + lOCs - 16 log^ 2 + 16n^ log^ 2 
+ ^7rMog2-P-lf^7r^-lp)C^C^ 

+ (3OC5 - SlTT^Cs - I5OC3 - fn' log 2 + ItT^ + if TT^ + f )C1 

+ (¥C3 + f vr^ W2 - f 7rMog2 - f vr^ + ^vr^ - 1^)C^T^ 
+ (4OC5 - Stt^Cs - I^Cs - f vr^ log^ 2 - If^vr^ log 2 

_i_ 28^4 I 4972 „2 3733 rp 

+ (-f|7r^ + ip)r^ + 4(6C3 + 7)C^i. 



+ 



(-67C3 + 607r2log2-f7r^-li)c| 

+ (-440(5 + IIIC3 + 347r2log2 + f tt^ + '-^)cICa 
+ (22OC5 - 2OIC3 - 327r' log2 + 6l7r' - ^)CfCI 
+ (220C5 + 175C3 + lf)C^l 



+ 



I /880 a I 1858\ ,^2t-i 1000 /-y rp2 3200 /^ t^2 

+ 264(C3 - 1)CffCa - 16(21C3 - 2)Cp^ 



1 

¥0 



+ 



4815 
16 



C^F + (-693(3 - 26471^ log 2 + IGStt^ + ^)c|C^ 



-(2 



+ (1205C3 - 367rMog2 - GOtt^ + ^)c|C 
+ (428C3 + 847r'log2 - IOStt" - ^)CfC\ + (-6I6C3 + '-f)C 
- 125CITf + (l327r2 - If^) C|C^T^ + (8471^ - ^)CfCIT^ 

4739 r^Srp I 4400 /-Y T^2 I 2800^2^^2 , f 3872 /- , 5324 \/-y ^2 
jg-O^iF i- ^^OpOA-t F + -^L^j^l F ' \ 3— S,3 -\ 9— j^FF^A 

I ( 4576 /- 1408 A /-I n j_ ( 704^ , 80\,^ 
+ (^~<53 9- j'^FA'^A + y 3-<,3 + -gj'^AA 

-^C% + 1656C|,Ca + (2178C3 + ^)c|C. 



r4 

F 



1 
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15 
'A 



(792C3 + ™)c|,Ci - (1386C3 + ^)C^Cl + '-fC 

1 



+ y^(l350C^ - 2664CfCA - 6140C|,C^ + 637Ci.C^ + 784C; 

- 3520CICaTf - 22A0CfCITf) {HCf + 7Ca) ^ 



+ ICf (-15C|, + 6CpC^ + 14C1) (llCp + 7C^)' § 

Po 



+ |CMllC^ + 7C^f ^, 

Po 



where 04 = Li4 , and 



^abcd^abcd ^abcd^abcd ^abcd^abcd 



(see notations in [43,45]). For SU{Nc) with Tp = i 



2' 



247VI ' 24 ' 24 



C Renormalization of dimension-4 operators 

As discussed in [31], three operators 

Ol = q (-iv ■ ^li)TK ^-iv-d {qflfVh^) 



(C.l) 



are closed under renormahzation for any F (note that 0\ — 0\ — id^ (qThy)). 
We have — Z{as{ijL))0^{ijL), where 



Z + 



Zb 




^OZaZ,^ 



, Z-' = Z-' + 



^„ Zb 
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7a 76 



7 = 7 + 



7a 7b 





If we take F = 7ar', then 

^10 ~ ^10 ' ^20 ~ ^20 ' ^30 ~ (3 ~ 2£)02o ' 



(C.2) 



(C.3) 



where 



(C.4) 



Of course, Oq = ■Z^(Q^s(/^))0'(/i). Wc obtain 



0[(/x) = 0;(/i) + Z (Z, + (3 - 25)Z,) 0^(/i) , 
02^(/x) = [l + Z (Z„ + (3 - 2£)Z6)] 0^(/x) , 
03r(//) = 30^(a^) . 



(C.5) 



Therefore, Z [Za + (3 — 2e)Zhj must be finite at £ — > 0. This aUows one to 
reconstruct 7^ from ^a' 

76 = -^(7a + A7a) , A7„ = i7„o(7ao-2/?o)(^)' + C»(«^). (C.6) 



The anomalous dimensions 'ja,b has been calculated in [31] with the two- loop 
accuracy (they are called 72,4 in this paper): 



a. 



+ 



10, 



-Tpni 



An 



+ 



The finite parts, at the next-to-leading order, are 



(C.7) 



olin)^o[in) + l^^o^oM, 



3 4n 



Ol{^) = 30M. 



1 M ^'^^^ 



o'M , 



(Ci 
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Note that 0[q = O'iq, but Of (/x) 7^ 0[{fi): additional counterterms in (C.5) 
yield a finite contribution, because of the 0{e) term in (C.3). 

As discussed in [32], two operators 

O'l^q {-iv • r/i„ , O^^iJ dxT{qrK, Ok(x)} 



are closed under renormalization for any F, and have 

Z = Z + 






i" 










+ 




:) 


, 7 = 7 + 




:) 















For F = 1, these operators are 02,3 (3.1); for F = 7°, they are Of ^ (4.1). This 
explains the relevant pieces of (3.3), (4.4). 

Three operators 



OT = -|g {-iv ■ <J^.F(1 + ^)a^''K , 
OT = -\q (-^^.) ^7M?^r(l + ^)a'^''K , 
O^^ij dxT{qTK,Om{x)} 



(C.9) 



are closed under renormalization for any F. Note that the indices /i, v live in 
the subspace orthogonal to v, due to -^hy — h^. These operators have [32] 



^ z 





Zb Z + Za 



ZZn 



7 = 7 + 



'^ ^ 

7b la 
\ lb la Im j 





Zb Z-^ + Za 



z-^z. 



\r7-\ 



(C.IO) 



The first two operators in (C.9) are (C.l) O3' and with F^ = — |i7''|^F(l + 
f)<7na'i therefore, the upper left 2x2 blocks in (C.IO) follow from (C.2). 

In the case F = 1, the operators (C.9) are related to (3.1): 

O- = -(l-^)(3-2£)O20, OZ = -{l-e)02o, OZ^O^. (C.ll) 
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We have 



+ [Z-'Z-'Z-^ - (1 - e)Z (Z- + (3 - 2e)Z^)] OM . 
Therefore, 

Z-'Z-'Z-^ - (1 - s)Z (Z^ + (3 - 2e)Z^ 



(C.12) 



must be finite at £ — > 0. This allows one to reconstruct 7™ in (3.3) from 7™^ 
in (C.IO): 



A7"^= 1 (7mo - 2/3o) (7^0 + 5%^) - haoTro ( ) + ^(«') 



47r 



The finite part, at the next-to-leading order, is 



(C.13) 



(C.14) 



In the case F = 7", the operators (C.9) are related to (4.1): 



= (1 - e)(l + 2e)0t, , = (1 - 2£)03"o + eOt, , 



We have 



(C.15) 



Otili) = Otiii) (C.16) 
+ [^-^Z-^Zr + (1 - 2e){Z + Z„)Z-] 03"(//) 

+ [^-^Z-^Z- + ((1 - 2£)Z, + eZ) + (1 - £)(1 + 2£)ZZr] 02{^^) . 
Therefore, 

z-x'^r + (i-2£)(^+^a)^r 



and 



Z-^Z-^Zr + ((1 - 2£)Z6 + eZ) + (!-£)(! + 2£)ZZ, 



m 
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must be finite at £ — > 0. This allows one to reconstruct 75^2 (4-4) from 7^^ 
in (C.IO): 



7r = 7r + A7r, 72"^ = 7r + A72", 

^iT = (7ao - 7mo + m (^)' + 0{al) , 



A7r = (7n.o - 2/?o) (7ro + 7.^) - i7ao7ro 



The finite part, at the next-to-leading order, is 



a 



47r 



47r 



(C.17) 



(C.18) 



The anomalous dimensions 71^2 has been calculated in [32] with the two-loop 
accuracy (they are called 7™]^^ in this paper): 



7r 



—2Cf- — \- Cp 



71 

9 



, 2 206\ 52^ 



etc 



a. 



47r 



9 3 ; 



10 2 46\ ^ 44^ 



etc 



The anomalous dimension 7"* (3.4) is, from (C.13) and (C.17), 



(C.19) 



r 



-7r - 372"^ + A7 , 



A7 = [(7mo - 2/3o) ilTo + 472^0) - |7ao7ro 



as 
An 



(C.20) 
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